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ABSTRACT. It is second attempt prediction results 
of classical biological control. We retrospective assesed 
the perspective of the introduction of two species of 
parasites (Blepharipoda scutellata R.-D. and Paraseti- 
gena silvestris R.-D.) against Lymantria dispar L. The 
difference between predicted level of steady-state den- 
sity of the host population and reality due to big errors 
estimated parameters (both for the parasites and the 
host). If for best investigated case the prediction is very 
approximate any assertion about absent perspective 
introduction parasites predators and competitor is false. 


РЕЗЮМЕ. Представлена попытка прогнозиро- 
вания результатов классического биологического 
контроля. Учитывая опыт прошлого, мы оценивали 
перспективу интродукции двух видов тахин 
(Blepharipoda scutellataR.-D. uParasetigena silvestris 
R.-D.) против Lymantria dispar L. Разница между 
прогнозируемым результатом уровня стабилизации 
(особая точка) популяции хозяина и действительно- 
стью, обусловлена большими погрешностями оце- 
ниваемых параметров (как для популяции паразита, 
так и для популяции хозяина). Несмотря на то, что в 
данном случае имеется относительно богатая ин- 
формация, результат все-таки получается очень при- 
близительным. Таким образом, всегда любое утвер- 
ждение относительно отсутствия перспектив интро- 
дукции энтомофагов является ненадежным. 


Introduction 


Since the first half of the 20-th century introductions 
of predator and parasitoid insects aimed at pest control 
were widely applied. The primery impulse was made by 
large-scale introduction projects carried out in the USA 
and Canada mainly to control gypsy moth and brown-tail 
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moth [Howard and Fiske, 1911; Mc Gugan and Coppel, 
1962]. Later on this methodology was used (with no 100 
% success) in controlling numerous agricultural and 
forest pests. Success has been reached in 16% and 
partially in 58% of the cases [Hall et al., 1980]. 

It should be mentioned that the methodology came to 
practice prior to be properly justified, since the theory 
had not been adequately developed. Perhaps, only a part 
of it considered the issues climatic suitability [Coppel 
and Mertins, 1977; Izhevsky,1990]was advanced. As for 
the biological basis of selecting parasitoid species for 
introduction, it has always been biased. In special liter- 
ature one can find a list of requirments that introduced 
species should meet [Coppel, Mertins, 1977; Hall et al., 
1980; Huffaker, Messenger, 1976; Izhevsky, 1990]. 
However, in practice low data quality or lack ofinforma- 
tion on some particular mechanisms of the introduction 
process did not slow down introductions. This approach 
was promoted, for instance, by Huffaker et al. [1971]. 
They wrote, in particular, that there is no need to study 
parasites and predators before introduction, since it is 
almost impossible to predict the effect of introduction. 
They proposed to introduce many parasitoid species so 
as to find out the best one. Probably, such an approach 
(supported to a certain extent by the majority of spe- 
cialists) along with the improperly planned and made 
experiments as well as lack of detailed observations on 
parasites and predators, provided relatively low suc- 
cess of the introduction projects. 

There is huge number of population dynamics mod- 
els for insects (literarute is not cited). However, any 
complete analysis of the introduction through even some 
promising models (see reviews Sharov, 1986; May, 
1976) has not yet been accomplished. No evaluation of 
the behavior of the models under conditions similar to 
the natural ones was done [Varley, 1974; Hassell, Var- 
ley, 1973]. We think that the main reason of this abnor- 
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mal situation is lack of publications on detailed biolog- 
ical observations (presented in numerical form) as well 
as low quality of some biological research. 

Computerised expert systems are powerful tool for the 
forecast of results of classical biological control. Expert 
system is simply an instrument. It is a good tool if it is 
comfortable to use and using it gives good results. Such 
systems are useful for the education and the practical 
biological control. Prognosis is based on the retrospection 
analysis experience classical biological control and the 
modern theory of dynamics of population density. 

But mathematical approach is typically not use in 
practical biological control programs. Perhaps because 
of many entomologists are lacking in higher level math- 
ematical knowledge. 

However, all the above does not mean that mathe- 
matical approach to the problem is impossible. Careful 
and critical review of published data as well as addition- 
al detailed observations on the majority of well-known 
pests and parasitoids can help provide valuable predic- 
tions using models. 


Methods 


Example of mathematical analysis of the introduction 

Just to illuminate an approach to prediction of the results 
ofintroduction we chose gypsy moth (Lymantria dispar L.) — 
one of the well-studied pests. Since detailed information on 
parasites and predators dynamics is needed to develop a model, 
we chose only two species parasitoids on which data were 
published in periodicals: 1) Blepharipoda scutellata R.-D. 
(=Sturmia scutellata R.-D., =Blepharipa scutellata R.-D., =B. 
pratensis Meig., =Blepharipoda pratensis Meig., =Herigia 
pratensisMeig.) (Diptera: Tachinidae); 2)Parasetigena sylves- 
tris R.-D. (=Phorocera silvestris R.-D., =Ph. agilis R.-D., 
=Ph. segregata Rond.) (Diptera: Tachinidae). Both species 
are solitary internal parasites of half-grown and larger gypsy 
moth caterpillars. Gypsy moth are the most preferable host for 
them [Zerova et al., 1989]. Both species have been introduced 
to the USA and well established there [Clausen, 1956]. 

Three data sets (7, 11 and 14 year time series) on the host 
and both parasitoids referred to the native area were extracted 
from publications [Znamenski, Lyamtcev, 1990; Panina, 1985; 
Lyamtcev, 1985; Sisojevic, 1975; Yafaeva, 1963]. These data 
were used as a basis for modeling. The authors mentioned 
above have used different units for the estimates of the 
population density: number per 100 cluster of leafs [Znamen- 
ski, Lyamtcev, 1990; Panina, 1985; Lyamtcev, 1985], number 
per tree [Lyamtcev, 1985; Yafaeva, 1963], number per 100 m? 
[Sisojevic, 1975]. All data were converted into number of 
individuals per 100 m?. The conversion was made using the 
methods proposed by Lyamtcev [1985] and regional tables of 
forest growth [Growth of the main ..., 1967]. 

Thus, we have only three rather short time series. Therefore, 
we need to simplify model structure as much as possible. The 
following assumption were used while constructing the model: 

1. Populations of host and parasites are local (neither 
emigration nor immigration are supposed). 

2. Both parasitoids act independently (1. e., probability of 
not getting attacked by both parasites is product of respective 
figures referred calculated for an each parasite separately). 

3. All hosts are accessible for both parasitoids (no 
shelters exist). 
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4. Parasitoids emerge from all the hosts attacked; mortality 
of parasitoids at the span “larvae in host — imago" is constant. 

Asa basis ofthe model we use logistic equation discribing 
host-parasite interactions by Nicholson and Bailey [1935] 
and Hassell and Varley [1969]. The model parameters were 
estimated using SYSTAT (version 4.1). 

The following equation shows how host population den- 


sity N,, at the next generation depends on its density N, at the 
current one: 

N,,,-N,*R*(1 - № /K) *EXP(-Q,,*P,, '»")* 

EXP(-Q,*P,,, n) (1), 


where R — fundamental net reproductive rate the host popula- 
tion, К — upper limit host population density, Р, — density 
of searching imago of B.scutellata of the current generation 
(number per 100 m?), P. a density of searching imago of P. 
silvestris of the current generation (number рег 100 т), О, — 
quest constant for B. scutellata, m, — mutual interference 
constant for В. scutellata, О — quest constant for P. silves- 
tris, m, — mutual interference constant for P. silvestris. 

It is also assumed that 
РІ = NO*EGG*0.5*WBL*S1 (2). 
P2 = NO*EGG*0.5*WPP*S2 (3), 
where NO*EGG*0.5 is number of half-grown gypsy moth 
caterpillars accessible for the parasites in previous genera- 
tion; 0.5 is survival of hosts at the stages egg-larvae of the 2- 
nd instar; № NS , Proportion of hosts infected by B. scutel- 
lata and P. silvestris in previous generation, respectively; P1 
and P2 post hunting density of populations of В. scutellata 
and P. silvestris, respectively. 

Lack ofaccurate data did not allow us to use variable EGG 
(mean number of eggs in the eggmasses in spring of the 
previous year). It was replaced by constant, namely by 350 
(mean number of eggs in eggmasses). 

The constants S1 and S2 were calculated as follows: 
51=0.45*0.5; S2=0.4*0.5, where 0.45 is survival of В. scutel- 
lata atthe period "larvae in host — imago" (the estimates given 
by Sisojevic [1975] and Zerova et al. [1989], were averaged); 
0.4 1$ survival of P. silvestris at the period "larvae in host — 
imago" (the estimates given by Sisojevic [1975] and Zerova 
et al. [1989], were averaged); 0.5 is proportion of females in 
parasitoid populations (actual data is not available). 

Dynamics of parasite populations was modeled as follows: 
pua ENS SI*U-EXP(-O. В, J) EXPO Е. m2) (4) 

Р м: N,S2*(I-EXP(- Q,*P.. Е —2))*ЕХР(- -Qp P Pi =!) (5), 
where Рик.» Pope — density of searching parasitoides of B. 
scutellata апаВ. silvestris of the next generaion, respectively. 

According to Hassell and Varley [1969]: 


Ор'" =Ln а (6) 

N,-W*N, i 
where N, — host population density at n-th generation, N — 
W*N, — density of hosts remained uninfected, W — propor- 
tion of hosts infected by both parasites. 

The estimation yielded the following values: 

B. scutellata: Q =0.213+0.022, m1=0.895+0.016, 
P. silvestris: Q,=0.048+0.032, m2=0.596+0.116. 

Net fundamental reproductive rate (R) = 1.28 

For the USA the data on the US populations was extracted 
from the article of Campbell [1967]. 

The role of the introduced parasites in 1910-1921 was 
considered as insignificant. 

The analysis of population dynamics of gypsy moth 
before introduction of parasites and predators was made using 
logistic model: N, ,,-N,*R*(1-N, /K) (7), 
where R= 1.580+0.408 and K= 347.704+66. 
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Results 


Prediction of the results of introduction was made on 
the basis of model (1). 

We assumed that the searching ability of parasites 
(expressed by Q and m) did not change during the 
introduction. 

For mean case, forecasted mortality causedB. scutel- 
lata and P. silvestris was 33% and 5%, respectively. 
Ferguson etal. [1994] reported mortality 32%, and 26%, 
respectively. Mortality due this species is 24% and 27%, 
respectivily in homeland. 

The model runs simulating introduction showed sus- 
tainable increase of parasitoid population densities after 
temporary gradual decrease. Usually they also showed 
that host population density declined while reaching 
stability (N*). Temporary increase of host population 
density may occur at the beginning of the process, 
especially if the magnitude of R is comparatively high. 
Oscillations of parasitoid population density occur only 
in some cases, for instance if m>1. The results of the 
simulation are presented at Table 1 and Fig.1. 

While running the model, the following ranges of 
parameters were investigated. As lower limits of Q, m 
(Fig.1) the minimum possible values (mean value minus 
error) were used; the only exception was made for О 
which lower limit was 0.001. The upper limits of the 
ranges were maximum possible values (mean value plus 
error). The ranges were devided into 10 (for Q and Q, m) 
and into 20 (for m and R). The runs were terminated 
when the difference between the current and subsequent 
values of the host density had become less than 0.0001. 
It usually took less than 50 iterations. 

The difference between levels of steady-state den- 
sity of the host population can be of four orders due to 
big errors in the estimated parameters (both for the 
parasites and the host). 


Discussion 


We know one attemt of prediction result of introduc- 
tion parasites [Varley, Gradvell, 1968]. 

The result of introduction is remarkable similar to 
the prediction for biological control in Canada made by 
Varley and Gradvell in 1968. 
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Fig.l. Probably steady-state density (№) of gypsy moth 
population past introduction parasites 

Рис. 1. Возможная точка динамического равновесия 
(устойчивая особая точка) популяции непарного шелкопряда 
после интродукции паразитоида. 


Flareups of winter moth in Nova Scotia and British 
Columbia have occured at about 9- to 10-year intervals 
(although not neccessarily on oaks, and not at the high 
levels predicted [Roland, Embree, 1995]. 

The steady-state density of gypsy moth population 
(№) is important characteristic of gypsy moth popula- 
tion dynamics. Most likely it is close to the mean 
population density. 

Analysis of the actual time-series showed that mean 
value of gypsy moth population density in the USA 
before application of classical biological control was 
about 58.0 eggmasses per 100 т. 

It dropped down to approximately 3.9 eggmasses per 
100 m? [Cambell, 1967] upon completing the first pro- 
grammes on introduction. If we also take into account 


Table 1. Predicted steady-state level of density of gypsy moth and its parasites population. 
Таблица 1. Предсказанная точка динамического равновесия для популяции непарного шелкопряда 


Companent of 
system 


12.4 


Bl. scutelata** 


10.5 





и его паразитоидов 


value of parametrs 


mean о, Q,, m, m, R Q,ts, Q,s, m+s, m24s, R 9-5 0,5, 5, 1—5 К о, Q,, m, m, В+ 


17.7 50.3 


* Mean number of eggmasses рег 100 т>, ** Mean number of pupariaper 100 т. 
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data by Bess [1961] the value droppes down more, to 
about 1.6 (logl10-scale)'. It is lower than threshold 
density that is recommended to switch on chemical 
control, 6 eggmasses per 100 n? , [Montgomery, 1990]. 
However, the problem of gypsy moth still remains par- 
tially due to inefficiency ofthe introduced parasites and 
predators [Price, 1990]. 

As it follows from our model analysis, introduction 
of two species of parasites (Blepharipoda scutellata and 
Parasetigena silvestris) may lead to stabilization of host 
population density. Its level (Fig. 1) is mainly higher 
than the threshold density recommended for switching 
on the control with chemicals. The predicted steady- 
state level of host population density is also much higher 
than the actual mean values. Partially И 1$ caused by 
artificially low amount of parasitoid specied (2) used in 
the model. Actually, there are about 10 established 
parasitoids and predators [Clausen, 1956] affecting gypsy 
moth. We also have not paid any attention to the mam- 
malian predators and birds. 

We have not made any comparison of the data on the 
role of the above two parasitoids after introduction with 
the predicted one. The reason 1$ that we were not able to 
find parallel long time series for the USA, characterizing 
Blepharipoda, Parasetigena and gypsy moth. There is no 
doubt that this data was collected many times, because 
many authors have published the analysis of these data 
[Burgess, Crossman, 1929; Bess, 1961; Campbell, 1967]. 

Theapplicability ofthe methodology described above 
will become finally clear by accomplishing the similar 
analysis with respect to other species. However, we 
think that is could be already accepted for preliminary 
evaluation of the results of introduction. 

Of course, the accurate data on host and parasite 
population dynamics in the native area are needed to 
succeed. 
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Appendix 1. Population density dynamics. 
Приложение 1. Динамика плотности популяции. 


Comments and reference 


Znamenski and 
Lyamtcev, 1990 
Lyamtcev, 1985 


Panina, 1985 


Yafaeva, 1963 








5.30 341.90 0.30 0.049 0.009 350.0* 

341.90 248.10 5.30 0.122 0.013 350.0* 

248.10 2.90 341.90 0.260 0.017 350.0* P. Sisojevic, 
2.90 0.06 248.10 0.397 0.386 350.0* 1975 
0.06 0.07 2.90 0.047 0.040 350.0* 





* Estimated values; density is given in number of individuals per 100 112; N1 — density of eggmass in spring of Nn-th year, N2 — density 
of egmass in spring of (Nn+1)-th year, NO — density of eggmass in spring of (Nn-1)-th year, Wbl — proportion of population parasitised 
by B. scutellata in (Nn-1)-th year, Wcc — proportion of population parasitised by P. silvestris in (Nn-1)-th year, Egg — number of eggs 
in eggmass in fall of (Nn—2)-th year or in spring of (Nn—1)-th year. 


